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The time evolution of the correlation functions of an ensemble of anharmonic N-component oscil- 
lators with 0(N) symmetry is described by a flow equation, exact up to corrections of order 1/N^ . 
We find effective irreversibility. Nevertheless, analytical and numerical investigation reveals that the 
system does not reach thermal equilibrium for large times, even when — > oo. Depending on the 
initial distribution, the dynamics is asymptotically stable or it exhibits growing modes which break 
the conditions for the validity of the expansion for large time. We investigate both classical 
and quantum systems, the latter being the limit of an 0{N) symmetric scalar quantum field theory 
in zero spatial dimensions. 
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Recently much progress was made in understanding 
the evolution of quantum field theories away from ther- 
mal equilibrium in the leading approximation [Q. 
However this approach has its shortcomings. The omis- 
sion of scattering leads to infinitely many conserved 
quantities [Q| which prevent thermalization. 

In we have taken the first step in an investigation of 
the evolution of an 0{N) symmetric scalar field theory 
in next to leading order The resulting system de- 

scribes the evolution of the 4-point couplings in addition 
to the 2-point functions of the fields and their conjugate 
momenta. At this level we can hope to answer the fun- 
damental question of whether the inclusion of scattering 
leads to thermalization. 

A full treatment of the non-linear behavior for the field 
theory requires much greater computational resources 
than at the leading order. One should therefore gain 
as much insight as possible in simple realizations of the 
problem In this letter we study the limit of zero 
spatial dimensions of an 0{N) symmetric (/j^-theory, or, 
equivalently the evolution of an ensemble of quantum an- 
harmonic oscillators. Our approach is based on the effec- 
tive action which generates the equal time one particle 
irreducible (IPI) correlation functions 

The time evolution of the effective action F obeys 
dtT[(p, TT;t] = — {Cci + £q) r[(p, tt; t], where the first op- 
erator, Cci, generates the classical evolution, whereas £q 
determines the dynamics of the quantum effects. We 
will deal with the specific example of a 0{N) symmet- 
ric ensemble of anharmonic oscillators Qa, a = 1, ...,7V, 
obeying the microscopic equations of motion dfQa — 
—rn?Qa — ^QbQbQa- The arguments of F are the mean 
values of the coordinates (pa = {Qa} and their conjugate 
momenta tto = (Qa) (in the presence of sources). From 
§,1), one infers (-0^ = (</)„, Tr^)) 
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On the level of IPI 2 and 4 point functions the most 
general form of F consistent with 0{N) symmetry is 



F = A^td), + B7r*TTa + 2C7r*0a} + ^[u(f>a(t)aMb 
+VTTa4>a4'b(l}b + WIT alT aMb + S [tT aTTb4'a<Pb ~ T^aT^aMb] 
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The time dependence of A{t), . . . , z{t) gives directly the 
time dependence of the propagator and the 4-point ver- 
tices. Higher vertices only induce corrections ~ 1/N'^. 
We emphasize that the truncation (^) still retains con- 
nected n-point functions of arbitrary order to the extent 
that they are constructed from 4-vertices and propaga- 
tors {Git = G5ab = {Qa{t)Qb{t)) - {Qa{t)){Qb{t)), etc.) . 

The exact equations for the inverse two point func- 
tions, with G = B/{AB - C^) and c = G/B, read 



dtG = -2cG, dtB = -2{c + -f)B, 
dtc = uj'^ - — ^ +JC + c?. 

tiLr 

Here, the time dependent frequency, uj, and 7 are: 
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For large N, with A, u, v, w, s, y, z scaling like l/A'^, one 
may expand in small XjN . We see that 7 ~ and lo 
becomes independent of u, w, . . . , in leading order. In the 
limit N ^ 00 the system of evolution equations for the 
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inverse 2-point functions is closed. In this approximation 
B/G becomes an additional conserved quantity |^. 

For a solution of at finite TV, however, one needs 
the time dependent 4-point functions. Their evolution 
involves, in turn the 6-point functions. The exact system 
is not closed and for any practical purpose we have to 
proceed to some approximation. Below we simply omit 
all contributions from IPI 6-point vertices. Then the 
evolution equations for the quartic couplings are 

dtu = uj'^v + iXBc - Xh^B^c^ - XBG'^\^2{N + 8)uc 
-{N + 8)vc^ + [{N + 2)w~ {N-l)s]c^Y 

df V = 2u?w - 4w - 7t; + 4AB - SXh^B^c^ - ABG^j 
2{N + 8)u + [{N - 10)w - 3(iV - l)s] + {N + 8)yc^Y 

dtw = 'iuj'^y - 3w - 27W - SXh^B^c - XBG^\^{N + 8)v 
+ [{N -10)w-3{N -l)s]c + 2{N + 8)zc^Y 

dts = 2uj^y - 2w - 27s - 2Xh^B^c - XBG^^{N + 6)v 

-2[4w +{N- l)s]c+ {N + 2)yc^ + Szc^j, (6) 

dty = ^uj"^ z - 2w - Sjy - Xh^B^ - XBG^ [{N + 2)w 

-(iV - l)s -{N + 8)yc + 2{N + 8)zc^^, 
dtz = -y - 'i'jz. 

Here we have omitted some terms ~ 1 /N'^. The neglected 
6-point functions can also be treated consistently as being 
~ 1/iV^. Thus, our set of evolution Eqs. (|^H) contains 
all contributions in the next to leading order The 
average energy e = ci + €2 = (H), 
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is conserved by the exact evolution equations, as well as 
by our truncation. Our system meets at least some basic 
conditions for thermalization: The number of interact- 
ing degrees of freedom N can be arbitrarily large and 
scattering effects are included. 

A first hint that the approach to equilibrium for large 
times is not obvious, however, comes from the study of 
possible time invariant solutions. The fixed points of the 
system Eqs. (|^^ correspond to stationary probability 
distributions. They obey, for uj^ > 0, 
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We observe a large manifold of fixed points since Eqs. (|^) 
have solutions for arbitrary B^ , s» and z^ ! This property 
seems not to be an artifact of the truncation. The fixed 
point manifold becomes even larger if one includes the 
six point functions and seems to characterize any higher 
(finite) polynomial truncation. Classical thermodynamic 
equilibrium, at temperature T (for h — 0) corresponds to 
the particular point in this manifold: 
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where lu'^^ obeys 
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with hmAT^oo t^oq = 5 ("^^ + Vrri^+2NXT) . 

The thermal fixed point in the quantum system occurs 
for non-zero In particular, we find in the quantum 

system a fixed point which differs from the classical equi- 
librium fixed point only by = A?i^i?^/(4w^), instead of 
^cq = 0. For this fixed point, all correlation functions in- 
volving only (f) are the same as for classical equilibrium. 
Inversely, for every fixed point of the quantum system 
there is one for the classical system, differing only by a 
shift in z*. Up to this shift, the quantum equilibrium 
could be reached by the classical system! 

An obstruction to thermalization arises from the ex- 
istence of an infinite set of conserved quantities. Be- 
sides {H), also arbitrary powers {HP{L'^Y) (with the 
squared generalized angular momentum) are conserved. 
For thermal equilibrium (L^) or (i?^) — {H)'^ are com- 
putable as functions of T, eg. (L^)^^ ^ N{N - l)rV^oq- 
Initial values of conserved quantities away from equilib- 
rium prevent thermalization. 

For classical statistics {Ti — 0), we have shown that 
the dynamical system describing the linearized motion 
orthogonal to the fixed point manifold is characterized 
by three doubly degenerate eigenvalues which are strictly 
imaginary. Consequently the motion is purely oscillatory. 
To probe the behavior in the non-linear regime, away 
from the fixed point manifold, we have solved the system 
of differential Eqs. (Qj^) numerically for various initial 
conditions. Unless otherwise stated, we consider below 
the the parameters N = 10, m? = 1, A = 1/N and ?i = 0. 

We first show, in Fig. |l|, the result of Gaussian initial 
conditions with G(0) = 0.18069, B(0) = 5 and aU other 
initial couplings zero. After a rapid initial increase of 
u the typical behavior becomes purely oscillatory. This 
confirms the results of our small fluctuation analysis and 
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the Fourier power spectrum reflects the relevant frequen- 
cies. On the other hand, the system never comes back 
to the initial Gaussian distribution which would require 
all quartic couplings to vanish simultaneously. Due to 
nonlinearities it clearly exhibits effective irreversibility! 
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FIG. 1. G{t) and u{t) for Gaussian initial conditions and 
h — 0. The second plot shows the Fourier power spectrum. 

The oscillatory behavior indicates that the thermal val- 
ues for the correlation functions can never be approached 
in a strict sense. Nevertheless, the time averages of 
some correlation functions over sufficiently large time in- 
tervals r could still approach the thermal equilibrium 
values. We therefore investigate quantities of the type 
{G)r{t) = {t)-^ Jl^ljldt'G{t'). Typical initial condi- 
tions for our study of the approach to equilibrium are 

G(0) =Geq(r)(l + A), ^/(0)=«eq(T), 

c(0) = u(0) = y(0) = w(0) = s(0) = z(0) = 0, (11) 



B(0)-i =T- AGeq(r) 
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Here T is the equilibrium temperature and GcqiT), 
Ueq{T) the corresponding values in Eqs. ( 



The ini- 
tial condition for B{0) is chosen such that the energy of 
the system e = ^ {3T + G^q{T)m^) , is independent of 
A. Thus, the initial conditions correspond all to the same 
effective temperature and the large time behavior should 
become independent of A if the system thermalizes. In 
Fig. 1^ we plot (G)ioo(i) for the initial conditions of (11). 
We see that for A not too large (G)ioo has for large t 
a stable mean value independent of r. The mean value 
depends, however, on the initial condition as specified by 



A. The curves correspond to the same effective T, but 
to different initial (L^). For A > we find a clear signal 
that the '2-point function' (Q^) never thermalizes. This 
constitutes the most important result of this letter. 

The situation for the quantum system is not very dif- 
ferent, as seen by comparing Figs. gA and ||C. For nega- 
tive = — 1 we show in Figs. ||b and ||D the results of 
the classical {h = 0) and quantum (?i = 1) evolution for 
T ~ 1 and with Gcq(r), Ucq{T) adapted appropriately. 
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FIG. 2. {G{t))ioo for initial conditions of (|ll|) with T = 1. 

We repeated the analysis of Fig. ^ for different N with 
similar results, even for very large N. The dependence 
of (G) on A becomes independent N as N ^ oo. 

The odd number of degrees of freedom and the absence 
of more conservation laws hints to the fact that the full 
nonlinear system of flow equations should not be Hamil- 
tonian, albeit being energy conserving. The spectrum 
of a Hamiltonian system is semi-positive definite and the 
system is stable in the sense that it has no growing modes. 
Otherwise this property is lost. Of course the existence of 
one growing mode implies others since the energy must 
be preserved by the evolution. This is precisely what 
happens to the full nonlinear system in the strong an- 
harmonic regime. As already visible in Fig. |2|thc system 
becomes unstable if the initial conditions arc far from the 
fixed point manifold. More precisely this means that the 
couplings grow large and the validity of the 1/iV expan- 
sion breaks down for large t. A useful criterion for the 
validity of the 1/N expansion requires that the contribu- 
tion to e proportional to u, |e„| = ^^^-^^^^^X\G^u\ remains 
smaller then ^ei. We denote by tc the time when this 
condition is violated and plot tc as a function of A in 
Fig.^for T = 1 and T — 5. The initial conditions are cho- 
sen here without the A terms in B{Q), i.e. B(0) = 1/T. 
For A < Ac the 1/A'' expansion remains valid at least 
for t < 10*. (The critical values are A^ = 0.655(1.012) 
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for T = 1(5) and ;i = and = 0.346(Ac < 10"^) for 
T — 1(5) and fi = \). For A > Ac one observes a very 
peculiar behavior of tc(A) with discontinuities at certain 
values of A. It suggests a chaotic dependence of the cor- 
relation functions on the initial conditions! (This type 
of behavior has been observed at leading 1/N j|]). This 
should not be confounded with the much more common 
chaotic behavior for individual solutions of the equations 
of motion. The latter, once averaged in an ensemble, of- 
ten leads to regular behavior of the correlation functions. 
Note that there are restrictions on the initial values of the 
correlation functions as they must be consistent with a 
positive probability distribution. We have not systemat- 
ically investigated these conditions so far. 
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FIG. 3. U vs. A, for T = 1,5 and B(0) = 1,0.2. The 
discontinuities seem to arise at arbitrary small steps in A. 
For A such that no point is plotted tc > 10^. 

Below we address the behavior for A > Ac in more 
detail in order to shed some light on the origin of the in- 
stability. Fig.^ shows a typical evolution for initial con- 
ditions generating instabilities. Notice that u{t) grows 
unbounded for large times while uj'^ (t) becomes negative. 
The power spectrum shows how the system goes stiff, 
particularly for the 4-point functions. In order to ensure 
the accuracy of our numerical integration we have used 
two independent methods, one based on an adaptative 
step algorithm and another especially tailored to handle 
stiff problems. Although the latter is intrinsically less 
precise we observed that instabilities occurred for coinci- 
dent tc- We conclude they are inherent to the dynamical 
system and robust to the choice of numerical integrator. 

From Fig. ^ it is apparent that many frequencies are 
shared among the 2 and 4-point functions. This may 
lead to secular instabilities since products of the 2-point 
functions act as sources for the dynamical equations of 
the 4-point functions. Secular instabilities result from an 
oscillator being driven at its own natural frequency and 
result initially in linear growth of its amplitude. For a 
given frequency uj shared by all degrees of freedom odd 
products of the 2-point functions can generate secular 



growth of the 4-point functions. In the quantum case 
these terms are more common, which may explain why 
the quantum motion generally becomes unstable sooner. 
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FIG. 4. The evolution of (^^{t), u{t) for initial conditions 
with T = 1, A = 0.666 and h — 1 and corresponding the 
frequency power spectra for G,c,u and z. 

In conclusion, we have presented a detailed study of 
the time evolution of correlation functions for ensembles 
of TV-component quantum and classical anharmonic oscil- 
lators. Our truncation method converges for large N and 
is well suited to describe an approach to thermal equilib- 
rium. Nevertheless, we find that the isolated systems do 
not thermalize! It will be very interesting to see if similar 
features persist in field theory or if effective thermaliza- 
tion occurs in the presence of many frequencies for the 
linear problem, due to dephasing . 
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